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Abstract 



The observed anisotropies of surface vibrations for unreconstructed FCC 
metal (110) surfaces are often reversed from the "common sense" expectation. 
The source of these reversals is investigated by performing ab initio density 
functional theory calculations to obtain the surface force constant tensors for 
Ag(llO), Cu(llO) and Al(llO). The most striking result is a large enhance- 
ment in the coupling between the first and third layers of the relaxed surface, 
which strongly reduces the amplitude of out-of-plane vibrations of atoms in 
the first layer. This also provides a simple explanation for the thermal contrac- 
tion of interlayer distances. Both the anisotropies and the thermal contraction 
arise primarily as a result of the bond topology, with all three (110) surfaces 
showing similar behavior. 
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The low-index surfaces of face-centered cubic (FCC) metals are arguably the most studied 
of surfaces. Though they are often intended to serve as prototypes for understanding the 
behavior of more complex surface systems, it has become evident in the last two decades that 
even these "simple" surfaces display a wide range of complicated and often counter-intuitive 
phenomena. In this paper, I consider the case of unreconstructed FCC (110) surfaces. I try 
to show that some of the observed features of their thermal behavior that have generally 
been accepted as violating "common sense" can actually be explained by simple arguments, 
backed up with results from ab initio density functional theory calculations. These results 
also have important implications for the study of other surface phenomena, and can, for 
example, provide insight into the mechanisms behind temperature-dependent surface phase 
transitions such as roughening and premelting. 

More bonds are broken to create a (110) surface from a bulk FCC crystal, than for a 
(100) or (111) surface, and one therefore expects the departures from bulk-like behavior 
to be largest for a (110) surface. Structurally, this is manifested in FCC metals in one 
of two ways: either the surface reconstructs into a "missing row" structure or the 
surface unit cell remains unchanged but there is a very pronounced oscillatory relaxation 
of interlayer spacings near the surface 0. In either case, the very existence of the surface 
lowers symmetry, and one expects to find anisotropic modifications in the thermal properties 
of atoms at or near the surface. 

The surprise is that for these unreconstructed FCC(llO) surfaces, some of these 
anisotropics are reversed, i.e., they violate the "common sense" expectation. For exam- 
ple, experiments show that in general: (i) atoms in the topmost surface layer have a bigger 
mean-squared displacement (MSD) in the surface plane than normal to it, whereas one might 
have expected the latter to be larger, since along this direction, there are no higher layers 
to repel the atoms in the first layer 0,^; (ii) the MSDs normal to the surface are larger 
for second layer atoms than for first layer atoms, though one would expect that the MSDs 
should decay monotonically into the bulk ||^ ; (iii) a third curious fact about unreconstructed 
FCC (110) surfaces is that in some cases there is a thermal contraction of interlayer distances 
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upon heating. The rule of thumb seems to be that if the first interlayer spacing di2 expands 
upon heating, then the next interlayer spacing ^23 contracts 0; however, if di2 exhibits a 
thermal contraction, then ^23 expands with increasing temperature P,^. 

Some of these patterns of behavior have also been observed in molecular dynamics (MD) 
simulations [0-0] but their origin remains a puzzle. In order to gain physical insight, and 
to be able to predict trends across differently oriented surfaces and different elements, one 
would like to not only reproduce this behavior in a theoretical calculation, but also to know 
where these anomalous properties arise from. Are they primarily a consequence of the bond 
topology, or are they due to peculiarities in the electronic structure? Such questions can 
reliably be answered by breaking up the process of obtaining a fully relaxed surface from 
the cleavage of a bulk crystal into steps, and performing a sequence of ab initio calculations 
at each component step. 

There have been few ab initio calculations of the lattice dynamics and thermal behavior 
of the (110) surfaces of FCC metals. Most of the previous calculations have been empirical 
(involving ad hoc modifications of near-surface interactions so as to reproduce experimen- 
tal phonon spectra) |]TU|-|T^, or semi-empirical The problem is, of course, that this 
procedure is not guaranteed to give a unique result , and one may be fooled into making 
modifications that are actually very far from the correct ones. The only FCC(llO) surface 
that has been investigated in detail, using first-principles methods, is the Al(llO) surface. 
The pioneering calculations of Ho and Bohnen [|1^ and other researchers ||15| , p!6| of the lat- 
tice dynamics of Al(llO) have recently been followed by an ab initio MD study 0, which 
reproduces the experimentally observed anisotropies and the thermal contraction of di2. But 
can the behavior of a simple metal like Al also serve as an adequate paradigm for transition 
or noble metals? Computational constraints make it currently unfeasible to perform such 
ab initio MD simulations on the noble metals or on transition metals. 

However, ab initio investigations of vibrational properties using the "frozen phonon" 
approach are possible, and have in fact been performed by previous authors on other low- 
index surfaces of these materials [jl^l- In this paper, I present the results of a series of such 



frozen phonon calculations on the bulk FCC structures as well as the unreconstructed (110) 
surfaces of Ag, Cu and Al. 

For each case, I have determined the exact interlayer and/or intralayer force constants 
from first principles, by distorting the equilibrium structure and computing the Hellman- 
Feynman forces |TB[ thereby induced on the atoms. These forces, as well as the relaxed 



structures, were obtained by performing density functional theory calculations using the 
package "fhi96md" |jl9| , with fully separable norm- conserving pseudopotentials and the 



Ceperley- Alder form of the local density approximation The calculations were carried 



out using plane wave basis sets with energy cut-offs of 60 Ry, 70 Ry and 20 Ry for Ag, Cu 
and Al respectively. It has been verified in earlier studies that the pseudopotentials and 
basis sets used here work well for both bulk and surface properties of these metals p2 



The surface calculations for Ag and Cu were carried out using a repeating slab geometry 
comprised of seven layers of atoms, separated by a vacuum thickness equivalent to five layers; 
while a 15-layer slab was used for the calculations on Al(llO). Bulk calculations were carried 
out using both the conventional cubic unit cell and a unit cell containing 12 to 16 layers of 
atoms stacked along the [110] direction. 

The interlayer force constants (pa/si'i j) were obtained by moving the atoms in a layer i 
along the direction a, and computing the force along (3 then experienced by atoms in other 
layers j. For the case of bulk layers stacked along [110], 1 find that for all three materials, 
the only significant elements of the interlayer force constant tensors are the diagonal terms 
(j)xx{h^^), (pyyihi ± 1) and (pzzihi ± 1) coupling nearest-neighbor layers, and (pz-Xhi ± 2) 
which couples next-nearest-neighbor layers (with the Cartesian directions defined as x = 
[110], y = [001] and z = [110]). All other elements of the interlayer force constant tensors 
are either zero by symmetry, or are smaller by at least an order of magnitude. Similarly, for 
the relaxed surfaces, the only significant terms involving the first (topmost) layer of atoms 
are ip^xil 2), (f)yy{l 2), ^^^(l 2) and 0^^(1 3). 

Some of these results are presented in Figure 1. For all three materials, I find that the 
elements of the interlayer force constant tensors stiffen considerably upon going from the 



bulk to the relaxed surface. The most notable feature is a dramatic increase in 02^(1 3), 
whose magnitude is almost doubled relative to the bulk value (pzzih'i =t 2). In fact, (pzzi^ 3) 
is found to be significantly larger than (pzzi^ 2). This means that if the atoms in layer 
1 are displaced along the z direction (normal to the surface) the resulting force along z 
experienced by atoms in layer 3 is considerably greater than the force on atoms in layer 2, 
which is a surprising and counter-intuitive result. Though the result 4>zz{^ 3) > 02^(1 2) was 
also obtained by some previous authors who fit force constant models to empirical data for 
Al(llO) [ill and Ni(llO) its importance seems to have been overlooked. I will argue 
that this large enhancement in 02z(l 3) is largely responsible for the anomalous thermal 
behavior of these surfaces. 

It should be emphasized that the results I have obtained for interlayer force constant 
tensors are exact, and do not involve any assumptions about the form or range of interatomic 
potentials. However, I now map these results onto a model potential, in order to better 
understand the implications of the changes in the surface force constant tensors. 

For each case, I have first assembled an extremely large database of results from ab initio 
frozen-phonon calculations. In addition to the terms listed above, the database includes 
results for other interlayer force constants, as well as intralayer terms. This database is then 
used to parametrize a simple form of interatomic potential: for each pair of nearest-neighbor 
(NN) atoms the interatomic potential U{r) is specified by a tangential parameter a = 
rQ'^[dU{r)/dr]r=ro and a radial parameter /? = [cPU (r) / dr'^]r=rQ, where rg is the equilibrium 
value of the interatomic distance r. 

For the interaction between two NN atoms in the bulk, I obtain {abb,f3bb} = 
{-.0007, .0181}, {-.0006, .0236} and {-.0007, .0152} for Ag, Cu and Al respectively, with 
all force constants being expressed in atomic units (Ha/bohr^). The oscillatory relaxation 
and charge redistribution at the relaxed surfaces result in a modification of these values. 
Accordingly, four kinds of NN bonds i-j between atoms in near-surface layers i and j (1-1, 
1-2, 1-3 and 2-3), are described by new parameters aij and (3ij, while all other terms are left 
unchanged from the bulk values. 
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Though this form of potential is admittedly simple, it is able to reproduce the gross details 
of the bulk phonon spectrum, and suffices to bring out the essential physics behind the altered 
surface behavior. Moreover, though the parameter set is heavily overdetermined (with the 
eight modified surface parameters being fit to a database of 34 different numbers determined 
from ab initio calculations for 13 different kinds of surface distortions, supplemented by three 
stability criteria), the quality of the fit is reasonable for Al, and extremely good for Ag and 
Cu. 

Next, the model potential is used to set up the dynamical matrix for a slab composed of 
many layers stacked along [110], which is then diagonalized to obtain phonon frequencies 
cUkA and eigenvectors ej^. The ath component of the MSDs at temperature T for atoms in 
layer i, {uj^{T)), is then given by p^ : 



where the summation runs over all wave-vectors k in the surface Brillouin zone and all 
phonon branches A; h is Planck's constant, M is the atomic mass and n^A is the Bose- 
Einstein distribution factor. 

In order to disentangle geometric effects from electronic ones, for each material I consider 
various cases. First, to determine the consequences of the reduced co-ordination at the 
surface alone, all NN interactions are replaced by the bulk parameters abf, and jSbf,. Upon 
using Equation (1), I find that all components of the MSDs of atoms in the first two layers 
are larger than the corresponding bulk values, but the anisotropics differ from those seen 
in experiments and MD simulations: the largest enhancement is in (uly) and {u1^), which 
are both approximately 2.2 times as large as in the bulk, and (ul^) > {u\^). Hitherto, there 
has been a tendency to attribute any anomalies in the behavior of FCC(llO) surfaces to the 
"very open" surface structure. However, these results show that the open structure alone is 
not sufficient to explain the observed phenomena. 

The situation is considerably altered for the three relaxed surfaces. The changes in the 
surface force constant tensors result in a considerable increase (15 to 50 %) in the value of 



the radial term Pu, and an even larger increase (45 to 85%) in /^is, relative to Pbb- The huge 
enhancement in /3i3 corresponds to the very large value obtained for 3), and implies 

that the bonds between NN atoms in layers 1 and 3 are extremely stiff. However, I find that 
the radial term coupling two NN surface atoms, Pu, is softened by ~20%. These results 
differ considerably from the 40% softening suggested for Pu for Ag(llO) by Franchini et al. 
13| , and the softening of [3i2 by 6% for Cu(llO) suggested by Black et al. who fit the 



parameters of their models to experimentally measured phonon spectra. 

What is responsible for the increased stiffness of the 1-3 NN bonds? To answer this, I 
looked also at the intermediate case of bulk-truncated surfaces (with electronic relaxation 
permitted, but all interlayer distances set equal to the bulk value). I find that though the 
value of /?i3 relative to /?fef, is slightly modified for the bulk-truncated surfaces, the huge 
enhancement comes upon going from the bulk-truncated surface to the fully relaxed one. 

The enhancement in /?i3 is sufficiently large to push {u\^) down significantly, now making 
{u\^) < {u\^) and {ul^) < {u\y). Thus, it is not that "the vacuum acts as a hard wall", as 
has been suggested as an explanation for the reversed anisotropics on Al(llO) 0, but the 
sub-surface atoms that constitute the hard wall that damps the vibrations along z of atoms 
in the surface layer. 

Another consequence of the increased value of /?i3 is that the surface will try to always 
maintain a fixed value for the interlayer separation (ii3, if necessary at the expense of changes 
in di2 and 6/23. This tendency has been confirmed by additional calculations in which, upon 
varying di2, ^23 was found to change in such a way that di^ was approximately constant. 
One can now understand why, upon heating an unreconstructed FCC(llO) surface, while 
di2 and ^23 may expand/contract, they usually do so in such a way that c/13 remains roughly 
constant, i.e., if one contracts the other expands. Of course a full treatment of the thermal 
expansion/contraction will require that one take into account the anharmonicity of the 
interatomic interactions; however, since the coefficient of thermal expansion is inversely 
proportional to the square of the harmonic force constant, the stiffening of /5i3 can be 
expected to have a big impact on the thermal variation of di^^. 



The fact that all three materials display the same trends, and the large impact of allowing 
for the relaxation of interlayer spacings, suggest that the enhanced stiffening of jSi^ over all 
other j3ijS may be more a consequence of the bonding geometry than of special features 
of the electronic structure. There are two relevant features in the geometry of FCC(llO) 
surfaces: (i) there is a very large reduction in the co-ordination of surface atoms, as a result 
of which the bulk-truncated surface relaxes by decreasing di2 significantly, (ii) a topological 
peculiarity of FCC(llO) surfaces [but not the (111) or (100) surfaces] is that a surface atom 
is connected by NN bonds to atoms in the first, second and third layers of atoms parallel to 
the surface. This explains why (pzzi^ 2) and ^22(1 3) may have comparable value, but not 
why the latter should be much larger. 

However, the large enhancement in /3i3 and the smaller increase in /?i2 can be rationalized, 
a posteriori, by simple trigonometry. One has to realize that for the FCC(llO) geometry, 
upon contraction of di2 (and thus also of dis), the shortening of the interatomic bond 
lengths and ri3 does not scale uniformly with the contraction of the corresponding 
interlayer separations di2 and di^. For example, a 10% contraction in di2, relative to the 
bulk interlayer spacing ds, translates to a 5% contraction in the NN bond length between 
atoms in layers 1 and 3, but a contraction of only 2.4% for the NN bond between atoms 
in layers 1 and 2. Since the stiffness of bonds scales as a high power of their equilibrium 
length [^, this results in a much larger increase in the radial force constant for 1-3 than 
1-2 bonds. 

The increase in ^23 that results from the oscillatory relaxation of these surfaces weakens 
this effect. However, these simple geometric considerations suggest that if ^23 is still small 
enough relative to ds so that c/23 + 2c/i2(i23 — 3c?^ < 0, one can still expect to find that 
/3i3 ^ /3i2. Similarly, if the contraction of ^34 is sufficiently small, then the radial force 
constant ^24 should be softened, thus increasing {u^z) further. 

Table I shows the results obtained for selected MSDs using the modified surface force 
constants and Equation 1. These results compare well with those deduced from measure- 
ments using low-energy electron diffraction (LEED), medium energy ion scattering (MEIS) 
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and helium atom scattering (HAS), as well as MD simulations using either ab initio (AI-MD) 
or Embedded Atom Method (EAM) potentials (there is, however, a considerable scatter in 
the values of MSDs available in the literature). It is important to note that my results 
do not include anharmonic effects, which may be large for FCC (110) surfaces, especially 
Cu(llO) Thus any discrepancy between my results and the experimental or MD ones 
may indicate that anharmonic effects are significant. 

The calculations presented in this paper can be extended by including anharmonic inter- 



actions within a quasiharmonic approximation, as has been done for other surfaces ||22| , |27 



However, I hope to have shown that the physics of the harmonic sector is itself interest- 
ing, and can go a long way towards predicting and understanding the thermal behavior of 
FCC(llO) surfaces. Given the simplicity of the geometric arguments, it should also be easy 
to extend these results to the unreconstructed (110) surfaces of other FCC metals. 
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FIGURES 

FIG. 1. Results from ab initio calculations for selected diagonal elements of the inter layer force 

constant tensors coupling tlie first layer of atoms witli sub-surface layers n, for the fully relaxed 
(110) surfaces of (a) Ag: xx,yy and zz elements (b) Ag, Al and Cu: zz elements. The zz elements 
for bulk Ag are also shown. 
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TABLES 

TABLE I. Comparison of calculated and experimental MSDs for atoms in the two topmost 

layers. T is the temperature in Kelvin, and (u^j) is the one-dimensional MSD of bulk atoms. 
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(Narasimhan, Fig. 1) 
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